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ABSTRACT 

We have observed two kinematically offset active galactic nuclei (AGN), whose ionised gas is 
at a different line-of-sight velocity to their host galaxies, with the SAMI integral field spectro¬ 
graph (IFS). One of the galaxies shows gas kinematics very different to the stellar kinematics, 
indicating a recent merger or accretion event. We demonstrate that the star formation associ¬ 
ated with this event was triggered within the last 100 Myr. The other galaxy shows simple disc 
rotation in both gas and stellar kinematics, aligned with each other, but in the central region 
has signatures of an outflow driven by the AGN. Other than the outflow, neither galaxy shows 
any discontinuity in the ionised gas kinematics at the galaxy’s centre. We conclude that in 
these two cases there is no direct evidence of the AGN being in a supermassive black hole 
binary system. Our study demonstrates that selecting kinematically offset AGN from single¬ 
fibre spectroscopy provides, by definition, samples of kinematically peculiar objects, but IFS 
or other data are required to determine their true nature. 
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1 INTRODUCTION 

In the event of a merger between two massive galaxies, the su¬ 
permassive black hole (SMBH) at the centre of each progenitor is 
expected to migrate towards the centre of the remnant (Begelman 
et al. 1980). When the two SMBHs become gravitationally bound 
to each other in a binary orbit, they continue to spiral inwards 
through the effect of three-body scattering of nearby stars (Quin¬ 
lan 1996; Khan et al. 2011, 2013) until at small radii (< 0.01 pc) 
gravitational wave emission causes further shrinking of the orbit, 
eventually leading to the SMBHs coalescing (Jaffe & Backer 2003). 

The remaining uncertainties in the mechanisms and timescales 
by which SMBHs coalesce, and their link to mergers which are a 
major driver of galaxy evolution, has led to a wide interest in iden¬ 
tifying and characterising SMBH binary systems. However, iden¬ 
tification of such systems has proved highly challenging. The only 
truly unambiguous cases are those in which two distinct sources 
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are observed in X-ray or radio interferometric imaging, indicating 
the presence of two accreting SMBHs (Owen et al. 1985; Komossa 
et al. 2003; Hudson et al. 2006; Rodriguez et al. 2006; Bianchi et al. 
2008). These systems will remain the ‘gold standard’ for SMBH bi¬ 
nary identification, but the difficulty of obtaining the necessary ob¬ 
servations, and the low redshift required to resolve the two sources, 
mean that significant sample sizes are currently beyond reach. 

Graham et al. (2015) and Liu et al. (2015) each claim detection 
of a SMBH binary system in a quasar via periodic modulation of the 
quasar’s brightness. The systems are intriguing and clearly warrant 
close investigation, but further confirmation is required before this 
method can be established as a viable technique for SMBH binary 
searches. 

A further technique for identifying SMBH binary candidates 
relies on the orbital motion of binary active galactic nuclei (AGN). 
If two SMBHs are in orbit with each other, and each carries a reser¬ 
voir of ionised gas travelling with the same velocity as the SMBH, 
then an optical spectrum of the system will, depending on orien¬ 
tation, show two sets of emission lines with different line-of-sight 
velocities, or double-peaked emission lines. Searches for such sys- 
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terns have produced many candidate SMBH binaries (Smith et al. 
2010; Liu et al. 2010), some of which have been confirmed by 
follow-up X-ray or radio imaging as described above (Comerford 
et al. 2011). Alternatively, long-term spectroscopic monitoring of 
candidates can identify those whose line-of-sight velocities change 
over time, indicating a binary orbit (Gaskell 1996; Eracleous et al. 
2012), although other variation in the spectra can mimic orbital 
velocity changes (Eracleous et al. 1997). Confirmation can also 
be found from optical/infrared imaging or, preferably, spatially re¬ 
solved spectroscopy, if two distinct sources are seen corresponding 
to the two kinematic components (Liu et al. 2010; Comerford et al. 
2011; McGurk et al. 2011; Fu et al. 2012; Villforth & Hamann 
2015). However, samples selected purely on the basis of having 
double-peaked emission lines are heavily contaminated by other 
systems that produce such spectra: Keplerian rotation, outflows, or 
other kinematic disturbances can all result in double-peaked lines 
that mimic SMBH binaries (Eracleous & Halpern 2003; Shen et al. 
2011; Fu et al. 2011). 

Recently, Comerford & Greene (2014, hereafter CGI4) pro¬ 
posed a related method to search for SMBH binaries in which only 
one SMBH is an AGN. In this case, a single set of emission lines is 
seen, but these lines are kinematically offset from the host galaxy 
redshift due to the same orbital motions that, in other cases, pro¬ 
duce double-peaked emission lines. CGI4 presented a sample of 
351 such objects drawn from the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) spectroscopic observations of galaxies. AGN with 
asymmetric emission lines, double-peaked emission lines, or mis¬ 
matched Baimer and forbidden lines, were not included in the sam¬ 
ple, in order to reduce contamination by outflows and other kine¬ 
matic features internal to the AGN. However, as noted by CGI4, the 
sample may still be contaminated by many of the same processes 
that can produce double-peaked emission, including outflows and 
dust obscuration. As a result, follow-up observations are required 
to confirm the candidates as SMBH binary systems. 

In this paper we examine in detail two of the kinematically 
offset AGN in the CGI4 sample, using integral field spectroscopy 
(IFS) from the Sydney-AAO Multi-object Integral field spectro¬ 
graph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an on¬ 
going survey of ^3400 low-redshift (z<0.12) galaxies with a wide 
range of masses and environments. For each galaxy, SAMI pro¬ 
duces spatially-resolved spectra over a 15" field of view, cover¬ 
ing the two wavelength ranges 3700-5700 A and 6300-7400 A at 
spectral resolutions of ~1700 and ~4500, respectively. The orig¬ 
inal SAMI instrument is described in Croom et al. (2012), recent 
upgrades and the sample selection for the survey are presented 
in Bryant et al. (2015), while Sharp et al. (2015) and Allen et al. 
(2015) give details of the data reduction and data products. 

The IFS data from the SAMI Galaxy Survey provide a useful 
test of the proposed binary SMBH nature of the selected galax¬ 
ies, primarily by supplying kinematic information about the stars 
and ionised gas throughout each galaxy. If the kinematic offset ob¬ 
served in the SDSS spectra is the result of the accreting SMBH 
being in a binary orbit with another, passive, SMBH, we would ex¬ 
pect the velocity of the AGN to be decoupled from that of the gas 
at larger radii. Away from the centre of the galaxy, the gas should 
in essence be unaware of the binary nature of the SMBH. As the 
galaxies here have been selected to have a kinematic offset at the 
centre, the key observational test is that the kinematics away from 
the centre should be normal. In other words, a binary system should 
result in a discontinuity in the observed gas velocity, with a spa¬ 
tially unresolved region at the centre of the galaxy in which there is 


a kinematic offset, surrounded by rotating gas that is settled in the 
same gravitational potential as the stars. 

Other potential causes of the kinematic offset would result in 
different spatially resolved signatures; some specific cases are de¬ 
scribed in Section 4. One important scenario is that of an outflow 
directed along the line of sight, which can produce the same sig¬ 
nature of a central kinematic offset surrounded by more normal 
gas that is expected from a binary SMBH. To distinguish between 
these two scenarios, a careful analysis of the properties of the cen¬ 
tral gas is required. Typical signatures of outflows include broad¬ 
ened and/or asymmetric emission lines, increased electron density, 
and shock heating. The presence of some or all of these features 
would indicate an outflow rather than a binary SMBH origin for 
the observed offset. Given that outflows are very common among 
luminous AGN, while SMBH binaries are expected to be rare, all 
possible signatures of outflows must be eliminated in order to gain 
a secure identification of a binary system. 

Throughout this paper we use a flat ACDM con¬ 
cordance cosmology: Hq = 70kms -1 Mpc -1 , = 0.3, 
Qa = 0.7. Unless otherwise specified, the emission¬ 
line transitions referred to are [Ol] = [Ol] A 6300; 
[On] = [O n] A3726 + [O n] A3729; [O ill] = [Om] A5007; 

[Nil] = [Nil] A6583; [Sn] = [Sn] A6716+ [Sn] A6731. Air 
wavelengths are used throughout. 


2 SAMPLE AND DATA 

The objects examined here are drawn from the CGI4 sample of 
kinematically offset Type 2 AGN. CGI4 selected these AGN from 
the OSSY catalogue (Oh et al. 2011) of emission and absorp¬ 
tion line properties, calculated from fits to galaxy spectra from the 
SDSS. The CGI4 objects were selected to lie above the Kewley 
et al. (2001) theoretical maximum starburst line in the [N II] BPT 
diagram (Baldwin et al. 1981), indicating the presence of an AGN, 
and to have a >3 <7 difference between the emission and absorp¬ 
tion line-of-sight velocities. AGN were rejected if their emission 
lines were asymmetric (weighted skewness ^0.5) as such lines are 
typical of outflows. 

Two of the CGI4 AGN have been observed as part of the 
SAMI Galaxy Survey, and in that context given the IDs 272831 
and 9016800393. These two galaxies were not targeted for their 
AGN properties, but rather they met the redshift and stellar mass 
cuts used for all galaxies in the survey, described in Bryant et al. 
(2015). Some basic information about the two galaxies, and details 
of the observations, are given in Table 1. The redshifts are measured 
from the SDSS spectra, while the stellar masses are calculated from 
the observed g and i magnitudes according to equation 3 in Bryant 
et al. (2015). Fig. 1 shows SDSS gri composite images of the two 
galaxies. 

The raw data were reduced using the standard SAMI data re¬ 
duction pipeline, described in Sharp et al. (2015). This pipeline 
carries out all steps from raw data to the production of fully cal¬ 
ibrated datacubes. The software itself is available online through 
its listing in the Astrophysics Source Code Library (Allen et al. 
2014); changeset 6c0c801, the same as for the SAMI Galaxy Sur¬ 
vey Early Data Release (Allen et al. 2015), was used to produce 
the data used in this paper. A full assessment of the quality of the 
pipeline’s results is given in Allen et al. (2015). 

The stellar kinematics of each galaxy were measured from 
the SAMI datacubes using pPXF (Cappellari & Emsellem 2004; 
Cappellari 2012), following the procedure described in Fogarty 
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Table 1 . Catalogue and observational information about the two kinematically offset AGN and the SAMI Galaxy Survey observations. 


SAMI ID 

272831 

9016800393 

SDSS name 

J120401.97+012641.5 

J011333.05+002948.0 

R.A. (J2000 deg) 

181.00821 

18.38774 

Dec. (J2000 deg) 

1.44487 

0.49667 

Redshift 

0.08352 

0.04404 

log 10 (Stellar mass / M 0 ) 

11.17 

10.63 

OSSY weighted offset (kms -1 ) 

67 ±13 

—74± 11 

Dates observed 

2013 Apr. 12, 13 

2013 Aug. 31; Sep. 3,4 

Number of exposures 

7 

7 

Total exposure time (s) 

12600 

12 600 

Seeing FWHM (") 

2.2 

2.1 


272831 9016800393 



-10 -5 0 5 10 -10 -5 0 5 10 


R.A. offset (arcsec) R.A. offset (arcsec) 

Figure 1 . SDSS gri composite images of the two galaxies examined in this 
paper, with IDs 272831 and 9016800393. The white circle in each panel 
shows the SAMI field of view. 

et al. (2014). Briefly, the 985 MILES stellar templates (Sanchez- 
Blazquez et al. 2006) were fit to a spectrum extracted from the cen¬ 
tre of the galaxy, to form a single master template. This master 
template was then fit to each spaxel individually, giving measure¬ 
ments of the velocity, velocity dispersion and higher moments at 
each location in the galaxy. Full details are given in Fogarty et al. 
(2014). 

In order to measure the emission-line properties from the 
SAMI spectra, the stellar continuum must first be subtracted. In 
principle this could be done using the fits derived during the stellar 
kinematics measurements, described above, but these are not opti¬ 
mised for accurate subtraction of stellar features in the spectra. To 
ensure that features such as the B aimer absorption lines are best 
removed, we re-fit the stellar continuum using a subset of 65 of 
the 300 available single stellar population (SSP) synthetic spectra 
from the MILES model library (Vazdekis et al. 2010), using the full 
subset at each spaxel. These model spectra form a regular grid of 
log(age) from 0.063 to 17.8 Gyrs and metallicity, [Z/H], from -1.31 
to 0.22 (equivalent to Z = 0.001 to 0.03). The 65 templates cover 
13 different ages at each of 5 metallicities; tests with the full set of 
50 ages showed no significant difference in the continuum fits. We 
include a multiplicative polynomial in the fit, in order to account 
for dust extinction and uncertainties in the flux calibration. Further 
details will be presented in Schaefer et al. (in prep.). 

Following subtraction of the stellar continuum, the emission 
lines were measured using FZIFU, an IDL package developed for 
the SAMI Galaxy Survey (Ho et al. 2014). LZIFU fits emission 
lines as the sum of one or more kinematic components, using the 
MPFIT implementation of the Levenberg-Marquardt least-square 
algorithm (Markwardt 2009). Each component consists of a Gaus¬ 
sian emission line for each of the requested transitions, with all 


transitions given the same velocity and velocity dispersion. Certain 
known flux ratios are also held fixed. To account for an uncertainty 
in the wavelength calibration of the blue arm of the spectrograph 
(Allen et al. 2015), an additional velocity offset is applied to all 
lines in that arm, with the magnitude of the offset left as a free pa¬ 
rameter in the fit. The result is a map of velocity and velocity dis¬ 
persion for each kinematic component, along with maps of the flux 
for each emission line in each component. The emission-line fluxes 
were corrected for Milky Way dust extinction using the Schlafly & 
Finkbeiner (2011) correction to the Schlegel et al. (1998) E(B — V) 
maps, and a Fitzpatrick (1999) extinction curve. 

For the centre of each galaxy we increased the signal-to-noise 
ratio (S/N) of the emission-line measurements by calculating the 
mean spectrum within a radius of 1", taking into account the corre¬ 
lations between nearby spaxels [see (Sharp et al. 2015) for a discus¬ 
sion of this issue]. The stellar continuum and emission lines were 
then fit again using the same procedure. 

We derived mass-weighted non-parametric star formation his¬ 
tories for our two objects using a spectral fitting technique follow¬ 
ing McDermid et al. (2015). Our approach will be described in de¬ 
tail in Scott et al. (in prep.), but we summarise it here. 

High S/N spectra were created by binning all spaxels within 
linearly-spaced elliptical annuli from the blue SAMI cubes for each 
object. Using pPXF, we fit to each spectrum 296 MILES SSP syn¬ 
thetic spectra with age from 0.0708 to 17.8 Gyr and [Z/H] from - 
1.71 to 0.22, again including a multiplicative polynomial. For each 
fit we imposed a smoothness constraint (using the REGUL key¬ 
word in pPXF) to the weights of neighbouring models. The relative 
weights assigned to each SSP model give the relative mass of each 
age and metallicity combination; we converted these to absolute 
masses by scaling each spaxel’s weights by the total stellar mass 
of the galaxy multiplied by the fraction of the galaxy’s continuum 
light (rest-frame 6700 A) in that spaxel. For the current work we 
summed over all metallicity values. The measurements presented 
here used only the blue arm of the SAMI spectra in the fit, but 
performing the same measurements on the combined blue and red 
spectra produced no significant difference in the results. 

Radio observations of the galaxies were obtained from the 
FIRST survey (Becker et al. 1995). We searched for radio sources 
within 30" of the coordinates of each galaxy, finding a match for 
272831 but not for 9016800393. For 272831 we confirmed the 
match by a visual inspection of an overlaid image. 

Both galaxies have FUV and NUV magnitudes available from 
the GALEX surveys (Martin et al. 2005; Morrissey et al. 2007). As 
272831 lies within the footprint of the Galaxy And Mass Assembly 
(GAMA; Driver et al. 2009, 2011) survey we obtained observed 
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Figure 2. Observed spectrum (red) and one-component LZIFU fit (black) for a central pixel from 272831. In the lower panels, seven wavelength regions 
around key emission and absorption lines are shown. These panels include the pPXF continuum fit (as used for the LZIFU emission-line fit) in green, and the 
wavelengths of key lines at the OSSY host galaxy (solid vertical lines) and emission line (dashed vertical lines) redshifts, as labelled on the plot. Note that the 
y range is different in each panel, and does not begin at zero. In all panels the grey dashed line shows the SDSS spectrum, scaled by the fibre area to place it 
on the same flux scale. 


magnitudes and k-corrections from the GAMA website 1 , and cor¬ 
rected them for Milky Way dust extinction using the same proce¬ 
dure as for the emission lines. For 9016800393 we corrected for 
Milky Way dust and then calculated rest-frame magnitudes using 
KCORRECT version v4_2 (Blanton & Roweis 2007). 

Neither galaxy is found in the ROSAT All-Sky Survey cata¬ 
logues (Voges et al. 1999), and other X-ray observations are not 
available. 


3 RESULTS 

As we will show in the following results, the observed properties 
and physical interpretations are very different for each of the two 
galaxies. Hence, we describe each one individually in the following 
subsections. 

3.1 272831 

The mean spectrum from the centre of 272831 is shown in Fig. 2. 
The lower panels show seven key wavelength regions in more de¬ 
tail, including both the observed spectrum and the pPXF/LZIFU 
fit. We use a single kinematic component to describe the observed 
emission lines, as a second component does not provide signifi¬ 
cant detections (S/N>3) of all of the strong lines. Vertical lines in 
the figure mark the positions of prominent emission and absorption 
features at the redshifts given in the OSSY catalogue. 

Fig. 2 also shows the SDSS spectrum for 272831. There is 
some disagreement between the overall shape of the SDSS and 
binned SAMI spectra, but this is expected given that the spectra 

1 http://www.gama-survey.org 


are formed from different apertures. The emission-line properties 
show a very close match between the two spectra. 

3.1.1 Velocity fields 

The emission-line and stellar velocity fields for 272831 are shown 
in Fig. 3, along with the difference between the two. The gas veloc¬ 
ity is taken from the LZIFU measurements, which simultaneously 
fit all strong emission lines with a single velocity. It is immediately 
apparent that the ionised gas and stars have very different kinematic 
properties: the stellar velocity field has very little rotation, but there 
is strong ordered rotation in the ionised gas. 

From the stellar properties we measure Ar, a parametrization 
of the specific angular momentum measured within one effective 
radius (Cappellari et al. 2007), as 0.100 ±0.014. Comparing to 
the Emsellem et al. (2011) analysis of the kinematics of early-type 
galaxies in the ATLAS 30 sample, 272831 is a slow rotator as mea¬ 
sured in the stars, lying below the cutoff of 0.121 for its observed 
ellipticity of 0.153. The values of Ar and ellipticity lie comfortably 
within the range observed by Emsellem et al. (2011) for slow ro¬ 
tators, which make up 14 per cent of all early-type galaxies. The 
ATLAS 30 sample shows a higher slow-rotator fraction for high- 
mass (above 10 11 M 0 ) galaxies like 272831; hence the stellar kine¬ 
matic properties of 272831, while uncommon in the context of the 
galaxy population as a whole, are not exceptional. 

The weak rotation of the stars is misaligned relative to the gas, 
with a difference in position angles of 66 ± 3°. The velocity field of 
the gas appears smooth, although we are not sensitive to structure 
on scales smaller than the seeing (FWHM 2'.'2). 

Fig. 4 shows the velocity as a function of position in the 
galaxy, taking a cut along the kinematic major axis of the ionised 
gas, approximately NE-SW. It is again clear that the gas and stars 
have very different kinematics - strong rotation in the ionised gas 
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Figure 3. Kinematic and emission-line flux maps for the galaxy 272831. From left to right, the top panels show: ionised gas velocity field; stellar velocity 
field; difference in velocities between gas and stars; log(Ha flux). The bottom panels show: log([N n]/Ha); log([S n]/Ha); log([0 i]/Ha); log([Om]/H/3). 
All velocities are in units of kms -1 , while the Ha flux is in units of 10“ 16 ergs -1 cm -2 pix -1 . Contours show the continuum flux at rest-frame 6400 A, as 
measured from the SAMI datacube, at intervals of 0.25 dex. In each of the first three panels, a dashed line shows the kinematic axis of the gas. In all panels the 
spatial coordinates are in arcseconds, with North up and East to the left. 



Figure 4. Velocity as a function of position in the galaxy for 272831. The 
two components plotted are stars (green circles, dashed line) and ionised gas 
(blue triangles, solid line). The vertical dashed line shows the centre of the 
galaxy. The cut is taken along the kinematic major axis of the ionised gas, 
with the positive direction lying approximately South-West. Dotted hori¬ 
zontal lines show the systemic velocity and, for illustration, ±300kms -1 . 


but not in the stars - but a further feature can now be seen. In the 
photometric centre of the galaxy, the ionised gas does not coin¬ 
cide in velocity with the stars, but is redshifted by 61 ± 8 km s -1 , 
in agreement with the OSSY measured offset of 67± 13kms -1 . 
Similarly, the velocities of the two flat regions in the ionised gas 
velocity curve, at ±325 ± 3kms -1 and —265 ±4kms -1 , are not 
symmetric around zero, but have a mean offset of 30 ± 3 km s -1 . 


3.1.2 Ionised gas properties 

Further panels in Fig. 3 show the distribution of Ha flux and key 
diagnostic line ratios across 272831. Fig. 5 shows the same data 
displayed in a series of ionisation diagnostic diagrams (IDDs), 
including also the velocity dispersion of the gas (Baldwin et al. 
1981; Veilleux & Osterbrock 1987; Rich et al. 2011). The central 
spaxels lie in the ‘LINER’ region of the [O lll]/H/3 vs. [N ll]/Ha, 
[O m]/H/3 vs. [S ll]/Ha and [O lll]/H/3 vs. [O l]/Ha diagrams, sug¬ 
gesting that the central gas may be ionised by shock heating (Kauff- 
mann et al. 2003; Kewley et al. 2006). Photoionisation by an AGN 
accretion disc, or shock heating in the immediate vicinity of such 
a disc, are disfavoured due to the spatial extent of the ionised gas: 
fitting a two-dimensional Gaussian to the Ha flux in the spaxels 
that lie above the Kewley et al. (2001) [N ll]/Ha classification line 
gives a FWHM of 3 '.'3, significantly larger than the 2 "2 seeing. 
Similar properties can also be caused by photoionisation by old 
stars (Binette et al. 1994; Singh et al. 2013), but this is unlikely 
to be the case in 272831 as the emission-line flux has a very dif¬ 
ferent distribution to that of the stars. In the outer regions of the 
galaxy the location on the IDDs is consistent with photoionisation 
by star-forming H II regions. 

Given the weakness of the H/3 emission line relative to the 
continuum level, it may be strongly affected by uncertainties in 
the template fit to the continuum. Such uncertainties would carry 
through to the [O m]/H/3 ratio and hence the IDD classification. 
However, we can place strong limits on the allowed range of H/3 
fluxes based on the strength of H/3 absorption observed in passive 
galaxies and the expected Balmer decrement (the Ha/H/3 ratio). 

The equivalent width (EW) in absorption of the H/3 line 
in the template fit to the central spectrum is 1.82 A, while the 
value measured of the observed spectrum (continuum and emis¬ 
sion line) is 0.58 A. If the template fit is inaccurate and the true 
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Figure 5. Ionised gas properties for 272831. Top row: line ratio diagrams for [N n], [S n] and [O I], from left to right respectively. Bottom row: line width as 
a function of [N n]/Ha, [S n]/Ha and [O i]/Ha, from left to right respectively. In each panel the size of the points is proportional to their Ha flux, and the 
colour shows the distance from the centre of the galaxy (red at centre, yellow at edge). In each panel, only points with S/N>3 in each of each of the lines 
featured in that panel are included. Cyan triangles on the line ratio diagrams show the OSSY measurements of the SDSS spectrum. The blue solid lines show 
the theoretical maximum starburst from Kewley et al. (2001), the green dot-dashed line shows the empirical starburst boundary from Kauffmann et al. (2003), 
and the red dashed lines show the Seyfert/LINER boundaries from Kauffmann et al. (2003, [N n] plot) and Kewley et al. (2006, [S n] and [O I] plots). The 
Kewley et al. (2006) classification for each region is labelled. The right-hand panel shows the mapping of Kewley et al. (2006) classifications back to the image 
of the galaxy: each spaxel is coloured according to its [N n] classification as Seyfert (dark blue), LINER (light blue), composite (light orange) or star forming 
(dark red). Contours are as in Fig. 3. 


EW is lower, the H/3 emission line flux will in turn be lower. 
However, the fitted value is already towards the lower end of 
the range observed in passive galaxies. At the most extreme, we 
can take the value of ~ 1.25 A observed in NGC4486 (Vazdekis 
et al. 2010), which would reduce the emission-line flux by a fac¬ 
tor of (1.82— 1.25)/(1.82 — 0.58) = 0.46, or an increase of 0.34 
in log([0 m]/H/3). Such a change would move the central points in 
Fig. 5 into the Seyfert region of the [O m]/H/3 vs. [N ll]/Ha dia¬ 
gram (top left panel), but would not change the classification in the 
other diagrams. 

An upper limit to the H/3 flux comes from the expected B aimer 
decrement of 2.86 (Osterbrock & Ferland 2006), for the extreme 
case of a galaxy with no dust reddening. Comparing to the observed 
value of 4.75, and noting that the Ha flux is much less sensitive to 
incorrect stellar templates, the true H/3 flux cannot be more than 
4.75/2.86 = 1.66 times greater than the measured value. Taking 
this limit would reduce the measured log([0 m]/H/3) by 0.22, in 
which case the central spaxels would still lie within the LINER re¬ 
gions in Fig. 5. In summary, despite the weakness of the H/3 emis¬ 
sion line we are able to place strong constraints on the line ratios, 
giving a confident classification of the galaxy. 

All the points in the three IDDs lie close to the fiducial shock 
line defined by Sharp & Bland-Hawthorn (2010), which traces the 
properties of a galaxy, NGC 1482, showing ionisation by shock 
heating. In NGC 1482 the shocks are associated with a starburst- 
driven outflow. There is an important difference in that the shock 
galaxies in Sharp & Bland-Hawthorn (2010) showed increasing 
[N ll]/Ha with increasing radius, while 272831 shows decreasing 


[N ll]/Ha with increasing radius, i.e. the spatial distribution is in¬ 
verted. This difference reflects the different sources of heating in 
the two galaxies: unlike 272831, NGC 1482 shows no signature of 
an AGN, indicating that the outflow is propelled and heated by star 
formation. Similar line ratios were also seen in a galaxy observed 
as part of the SAMI commissioning (Fogarty et al. 2012), and in¬ 
terpreted as a shock heating due to an outflow driven by star for¬ 
mation; again the spatial distribution of ionisation properties was 
reversed relative to 272831. 

A further galaxy from the SAMI Galaxy Survey with related 
ionisation properties is 209807 (Ho et al. 2014), again interpreted 
as star-formation-driven shock heating. This interpretation was re¬ 
inforced by a strong correlation between the ionised gas velocity 
dispersion and the [N ll]/Ha, [S ll]/Ha and [O l]/Ha ratios, which 
is also seen in 272831, as shown in the fourth, fifth and sixth panels 
of Fig. 5. In 209807 the higher ionisation gas was detected in the 
centre of the galaxy, but only in the form of a distinct kinematic 
component forming the outflow. We do not detect a separate com¬ 
ponent of this sort in 272831; a single component is sufficient to de¬ 
scribe the emission across the entire galaxy. In summary, although 
the IDDs are consistent with shock-heated gas, the kinematic in¬ 
formation suggests a different configuration to the outflows seen in 
other similar galaxies. 

We measure the gas-phase metallicity of 272831 us¬ 
ing the combined [Nll]/Ha and R 23 = ([Oll]A3727 + 
[O III] AA4959,5007)/H/3 calibration of Mannucci et al. (2010), in 
turn based on the calibrations of Maiolino et al. (2008). We choose 
this calibration in order to allow direct comparison with the Man- 
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Figure 6. Upper panel: SFR as a function of age, as measured from the 
stellar populations, for the entirety of 272831 (solid blue line) and only the 
outer regions (dashed green line). Lower panel: Cumulative mass fraction 
as a function of age, i.e. the fraction of stellar mass in stars younger than 
the x-axis value. 

nucci et al. (2010) mass-metallicity relation. In the outer regions 
of 272831, where the emission-line properties are consistent with 
star formation, the metallicity is 12 + log(0/H) = 8.93 ± 0.04, 
with no obvious spatial variation. Using each of the [N ll]/Ha and 
R 23 calibrators in isolation gives similar values: 8.87 ±0.03 and 
8.98 ±0.02, respectively. 

We estimate the Eddington ratio of the AGN from the ob¬ 
served [Olll] luminosity, using the [O lll]-to-bolometric conver¬ 
sion L\y oi = 3500 L [ 0 hi] from Heckman et al. (2004). Extracting the 
[O III] luminosity from a circular aperture with radius 2" centred 
on 272831, we find 0 i = 1.7 x 1O 1O L 0 , with an uncertainty of 
0.38 dex due to the scatter in the Lboi _ ^[om] relation. Note that this 
luminosity, and hence the derived Eddington ratio, is likely an over¬ 
estimate if some or all of the central [O III] flux arises from shocks 
or other ionisation sources. We measure a central stellar velocity 
dispersion, <j v , of 249 ± 6 km s _1 , and use the Giiltekin et al. (2009) 
relation to convert this to a black hole mass, Mbh> of 4.0 x 10 8 M 0 . 
The scatter in the M^\\-o v relation gives an uncertainty of 0.31 dex 
in Mbh- Combining the Lbol and Mbh measurements, we derive an 
Eddington ratio of 1.3 x 10 -3 , with an uncertainty of 0.49 dex. 


3.1.3 Stellar populations 

The upper panel of Fig. 6 shows the star formation rate (SFR) as a 
function of age for 272831, measured from the stellar populations 
fits. The lower panel shows the cumulative stellar mass as a fraction 


of the total, plotted as a function of age. The vast majority of the 
stellar mass was formed more than 1.6 Gyr (10 92 yr) ago, with a 
small burst of star formation occurring within the last 0.1 Gyr. This 
recent burst has formed less than 0.1 per cent of the total mass of 
the galaxy. 

In order to compare the stellar populations results to the cur¬ 
rent SFR measured from Ha we discard the central region, up to a 
radius of I'.'l, and retain only the outer regions of the galaxy. Do¬ 
ing so avoids contamination by the AGN or shocks. The discarded 
region is 1.5 times the size of the seeing disc for the data, and cor¬ 
responds to a physical size of 2.7 kpc. It encloses the region of high 
[N ll]/Ha seen in Fig. 3 so can be assumed to contain the majority 
of emission from gas ionised by sources other than star formation. 
The remaining emission in the retained region is dominated by star 
formation. After correcting for dust extinction using the observed 
Baimer decrement and the jUo .4 dust curve from Wild et al. (2011), 
we sum the extinction-corrected Ha flux in this region. We then 
convert this flux to an SFR using the relation in Kennicutt (1998), 
converted to a Chabrier (2003) initial mass function (IMF), giv¬ 
ing 0.66 ± 0.04 M 0 yr -1 . If the Ha flux across the entire field of 
view is assumed to be due to star formation, the SFR increases 
to 1.05 ± 0.05 M 0 yr _l , but this is certainly an overestimate due 
to contamination by other ionisation sources in the centre of the 
galaxy. 

The Ha SFR is close to the SFR inferred from the youngest 
stars in the stellar populations fit: a mass of 4.0 x 10 7 M 0 was 
found for the population with age 71 Myr, corresponding to an SFR 
of 0.56 M 0 yr -1 averaged over the age of the bin. In this case no 
young stellar populations were found in the discarded central re¬ 
gion of the galaxy, so the SFR is the same for the outer regions as 
for the entire galaxy. 


3.1.4 Radio emission 

The galaxy is detected in the FIRST radio catalogue, with an in¬ 
tegrated flux at 1.4 GHz of 4.27 mJy, which corresponds to a radio 
luminosity of 7.54 x 10 22 WHz -1 . Using the calibration of Hop¬ 
kins et al. (2003) and assuming that all the radio continuum emis¬ 
sion is due to star formation, this would correspond to an SFR 
of 41.7 M© yr -1 , far higher than the upper limit of 1.05 M 0 yr -1 
implied by the Ha emission. The vast majority of the radio con¬ 
tinuum must therefore be due to an AGN. The emission lines in 
the centre of the galaxy have an excitation index, as defined by 
Buttiglione et al. (2010), of 0.55 ±0.02, below the cutoff of 0.95, 
making 272831 a low excitation radio galaxy (LERG). This clas¬ 
sification is consistent with the radio luminosity, which is in the 
range dominated by LERGs (Best & Heckman 2012), and with the 
low Eddington ratio of the AGN. 


3.1.5 UV emission 

The GALEX observed magnitudes for 272831 are 20.51 ±0.09 and 
19.65 ±0.04 for the FUV and NUV filters, respectively. After re¬ 
moving the effect of dust extinction within the Milky Way, and k- 
correcting, we derived rest-frame magnitudes of 20.23 ± 0.09 and 
19.12 ±0.04. Using the calibration of Salim et al. (2007) for ‘nor¬ 
mal’ galaxies, again assuming a Chabrier (2003) IMF, these corre¬ 
spond to an SFR of 12.8 ± 3.0M 0 yr -1 . This calibration includes 
an empirical correction for intrinsic dust extinction, based on the 
FUV—NUV colour. The quoted uncertainty does not include the 
uncertainty in this correction, which increases the flux by a factor of 
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Figure 7. As Fig. 2 but for the galaxy 9016800393. In this case a two-component LZIFU fit is used; each emission-line component, added individually to the 
stellar component, is plotted in blue. The lower panels include the LZIFU fits to each component. 


37, so will be an underestimate. However, the large dust correction 
is consistent with the relatively high Balmer decrement, 4.8 ± 0.2, 
seen in the centre of the galaxy. 

The SFR calculated here may be contaminated by UY flux 
from the AGN, but Salim et al. (2007) show that such contami¬ 
nation should not exceed 15 per cent of the total flux, based on 
an extrapolation of the limit on r-band contamination measured by 
Kauffmann et al. (2003). Conservatively assuming the maximum 
level of contamination gives a final UV SFR of 10.9 ± 2.6 M 0 yr _1 . 
This value is an order of magnitude higher than the Ha SFR, which 
may be due to the different timescales probed by the two methods. 
The Ha flux is measuring star formation more recent than the peak 
seen in Fig. 6, but the UY flux, which probes timescales up to a few 
10 8 yr, includes that peak and so is expected to give a higher SFR. 
The discrepancy between the two measures may also be enhanced 
by uncertainties in the dust correction, which can vary greatly be¬ 
tween individual galaxies (Boselli et al. 2009). 


3.2 9016800393 

The central spectrum of 9016800393 is shown in Fig. 7. In this 
case, two velocity components are required to provide a good fit 
to the emission lines: a narrow component (<T V = 107 ± 1 kms -1 ) 
with a blueshift of 44 ± 2 km s - Relative to the stars, and a broader 
component (<r v = 304 ± 7 km s -1 ) that is blueshifted by a further 
229 ± 14kms -1 relative to the narrow component. The broader 
component is measured at high S/N, formally an 18-<7 detec¬ 
tion in the [O III] A5007 line. The blueshift measured by OSSY, 
74 ± 11 km s -1 , lies between those for the two components in the 
SAMI data, suggesting it was influenced by both physical compo¬ 
nents. 


3.2.1 Velocity fields 

The stellar and ionised gas velocity fields for 9016800393 are 
shown in the upper panels of Fig. 8. In the ionised gas we show 
for simplicity a single-component fit to the data. Even in the cen¬ 
tral region, the velocity measured in this way is dominated by the 
narrower of the two velocity components. 

In contrast to 272831, both the ionised gas and the stars in 
9016800393 show simple ordered rotation, with position angles 
that agree within the measurement uncertainties. In the outer re¬ 
gions of the galaxy (radius 5-7"), the velocity of the stars has a 
mean value of 79 per cent of that of the ionised gas, using the 
single-component LZIFU fits. This offset is as expected due to 
asymmetric drift, the tendency of stars to lag behind the gas in ro¬ 
tation velocity due to their increased random motions. The value 
observed in 9016800393 is close to the mean value of 75 per cent in 
the SAMI Galaxy Survey sample as a whole (Cortese et al. 2014), 
and comparable to the 89 ± 8 per cent measured by Martinsson et al. 
(2013) from 30 spiral galaxies in the DiskMass survey. 

Fig. 9 shows the velocity as a function of position in the galaxy 
for each of the components, taking a cut along the kinematic ma¬ 
jor axis. Each of the two components of the ionised gas is shown. 
These two components have markedly different behaviour: the nar¬ 
row component traces the stellar kinematics closely, other than 
the asymmetric drift described above, but the broader component 
shows a blueshift of ~250kms -1 in the centre of the galaxy, as 
seen in the single spectrum shown in Fig. 7. 

A careful examination of the spectra in Fig. 7 shows that the 
narrow component of the [O III] emission lines is blueshifted rela¬ 
tive to the same component in other emission lines, with an off¬ 
set of 71kms -1 , bringing its total blueshift relative to the stars 
to 115kms -1 . The LZIFU fitting procedure has incorporated this 
shift within the parameter that defines a wavelength calibration off¬ 
set between the blue and red arms of the spectrograph. However, 
in this case the excellent agreement between the SAMI and SDSS 
emission line positions in both arms indicates that the offset is not 
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Figure 8. As Fig. 3, for the galaxy 9016800393. The measurements from the single-component LZIFU fit are shown. 



Figure 9. Velocity as a function of position in the galaxy for 9016800393. 
The components plotted are: stars (green circles, dashed line), narrow gas 
(blue triangles, solid line) and broad gas (black stars, solid line). The verti¬ 
cal dashed line shows the centre of the galaxy, while the horizontal dotted 
line marks the systemic velocity. The cut is taken along the kinematic major 
axis, with the positive direction lying approximately East. 

due to a calibration error but a genuine shift in the [O III] velocity. 
The offset is not apparent in the other forbidden lines. 

3.2.2 Ionised gas properties 

Maps of the emission-line properties in 9016800393 are shown in 
Fig. 8, using the single-component fits to the emission lines. The 
upper panels of Fig. 10 show the IDDs for the same fits. In the lower 
panels, the same diagrams are shown with the emission separated 
into the two kinematic components observed. 

Looking first at the single-component fits, it is clear that the 
emission lines in the centre of 9016800393 are dominated by a 


Seyfert-like AGN, with some star formation in the outer disc. The 
points lie close to the fiducial AGN-excited line defined in Sharp 
& Bland-Hawthorn (2010). The l-<7 width of the lines does not 
exceed 150kms -1 across the galaxy. 

When two components are fit to the emission lines, as shown 
in the lower panels of Fig. 10, the narrower component shows es¬ 
sentially the same results as for the single-component fits. The 
broad component in the outer spaxels also follows the single¬ 
component results, suggesting that in this region it may be due to 
seeing effects and integrating along the line-of-sight in an edge-on 
galaxy producing non-Gaussian emission line profiles, rather than 
being a separate physical component. However, in the central re¬ 
gion the broader component shows some marked differences: it is 
significantly broader than any of the single- or narrow-component 
fits (~300kms -1 as opposed to <150kms _1 ), the [Om]/H/3 ra¬ 
tio is somewhat higher, and it has a lower [S ll]/Ha ratio (loga¬ 
rithm —0.6, compared to —0.2 for the narrow component in the 
same region, lower by a factor of ~4). The broader component is 
marginally spatially resolved: a 2D Gaussian fit to its spatial pro¬ 
file in [O III] gives a FWHM in the major/minor axes of 3"4/2. /r 6, 
compared to a seeing FWHM of 2'.'I (l" corresponds to 876 pc at 
the redshift of 9016800393). The extension of the broader compo¬ 
nent is aligned 11° anticlockwise of the E-W direction, or 18° away 
from the major axis of the galaxy. 

The unusually low [Sll]/Ha ratio for material excited by 
an AGN suggests the broad component may be collisionally de- 
excited, indicating that the electron density is at or above the critical 
density for the relevant transitions, 1.2 x 10 3 and 3.3 x 10 3 cm -3 
for the 6716-A and 6731 -A transitions, respectively (Rubin 1989). 
However, the limited S/N of the broad component measurement, 
which will be compounded by any uncertainties in the continuum 
subtraction, prevent us from placing strong constraints on the den¬ 
sity. 

The log([N ll]/Ha) ratio in the outer regions of 9016800393 
is —0.29, higher than the upper limit for reliable metallicity mea¬ 
surements in gas ionised by young stars (Mannucci et al. 2010), 
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Figure 10. Ionised gas properties for 9016800393. The upper two rows are as in Fig. 5, showing the IDDs for single-component fits to the emission lines. The 
panel on the right shows the mapping of [N n]-based classifications back to the galaxy, using the same colour coding as in Fig. 5. The lower two rows show 
the results for the two-component fits: the red circles show the narrow component and the blue triangles show the broad component. In each panel, each line 
in the diagram must have S/N>3 in the individual component for a point to be included. 


so we use only the R 23 calibration. The resulting metallicity is 
12 + log(0/H) = 8.97 ±0.07. 


Following the same procedure as in Section 3.1.2 to derive 
the Eddington ratio for the AGN in 9016800393, we measure 
Lboi = 8.4 x 10 9 L 0 and Mbh = 1.5 x 10 7 M 0 . These combine to 
give an Eddington ratio of 1.7 x 10 -2 , with an uncertainty of 
0.58 dex due to the scatter in the scaling relations. 


3.2.3 Stellar populations 

The SFR and cumulative mass fraction for 9016800393, from the 
stellar populations fits, are shown in Fig. 11. The shape of the SFR 
distribution as a function of time is similar to that of 272831, with 
an extended period of star formation at ages above 1.6 Gyr, and a 
strong peak of recent star formation within the last 0.13 Gyr. The 
recent star formation accounts for 0.3 per cent of the total mass of 
the galaxy. 

Again discarding the central region of the galaxy, we find a 
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Figure 11. As Fig. 6, for the galaxy 9016800393. 


mass in the youngest stars (71 Myr) of 4.3 x 10 7 M 0 , correspond¬ 
ing to an averaged SFR of 0.61 M 0 yr -1 . The Ha SFR for the 
same region is substantially lower, 0.060 ± 0.007 M 0 yr -1 . Inte¬ 
grating across the entire galaxy gives a mass of 5.5 x 1O 7 M 0 , 
equivalent to an averaged SFR of 0.78 M 0 yr -1 , and an H a SFR 
of 0.263 ± 0.019 M 0 yr -1 . However, using the entire galaxy will 
strongly contaminate the Ha SFR with the AGN emission. 

3.2.4 Radio emission 

There is no detection in FIRST for 9016800393, with a cata¬ 
logue detection limit (including clean bias) at its coordinates of 
0.86 mJy, putting an upper limit on the radio power of the AGN 
at 4.0 x 10 21 WHz -1 . This limit is much higher than the radio 
power expected due to star formation, which lies in the range 1.1- 
4.8 x 10 20 WHz -1 , depending on the fraction of Ha emission aris¬ 
ing from the AGN. 

3.2.5 UV emission 

The GALEX observed magnitudes for 9016800393 are 21.56 ± 
0.08 and 20.80 ± 0.07 for the FUV and NUV filters, respec¬ 
tively. After removing the effect of dust extinction within the 
Milky Way, and k-correcting, we derived rest-frame magnitudes 
of 21.32 ±0.08 and 20.48 ±0.07. Using the calibration of Salim 
et al. (2007) for ‘normal’ galaxies, these correspond to an SFR of 
O.54±O.15M 0 yr -1 . As before, the quoted uncertainty does not 
include the uncertainty in the intrinsic dust correction, so will be 
an underestimate. Correcting for contamination by the AGN gives 
a final UV SFR of 0.46 ± 0.13 M 0 yr“*. 


4 DISCUSSION 

Again, due to the very different properties of the two galaxies re¬ 
vealed by the SAMI observations, we examine each one individu¬ 
ally in the following subsections. 


4.1 272831 

The striking difference between the gas and stellar kinematics in 
272831 implies an external origin for the gas, either from a recent 
merger with a gas-rich galaxy, or direct accretion. In either case the 
stellar populations indicate that the star formation associated with 
the event occurred within the last 100 Myr. This age may not corre¬ 
spond directly to the timescale for the event itself as, for example, a 
merger may trigger star formation at first passage, as much as 1 Gyr 
before the final coalescence (Patton et al. 2013). The UV SFR is 
higher than either the Ha SFR or the SFR inferred from the stel¬ 
lar populations, which may result from the uncertainty in the UV 
dust extinction correction, or simply the scatter in the relationship 
between UV luminosity and SFR. 

Comparing to fig. 1 of Mannucci et al. (2010), the metallicity 
of the gas, 12±log(0/H) = 8.93 ±0.04, is somewhat low for a 
galaxy with a stellar mass of 10 1117 M 0 , lying at the edge of the 
90 per cent range illustrated in the figure. The median metallicity 
at that mass is 9.16. The observed metallicity is more typical of 
a galaxy with a stellar mass in the range 10 9 7 -10 10 - 2 M 0 , possi¬ 
bly implying that the gas results from a merger with a galaxy of 
approximately this mass. 

The kinematic centre of the gas is offset in velocity from the 
stars, resulting in the offset detected by CGI4. This offset implies 
that the gas has yet to settle into a stable configuration, in agreement 
with the relatively short timescale for the accretion or merger event. 
The ionisation properties in the centre are consistent with shock 
heating, indicating that the ionisation may not be directly due to 
photoionisation from an AGN. Instead, it may be an indirect result, 
due to interaction of an AGN jet with the gas, or not related at all. In 
this last possibility, the gas may be shock heated due to interaction 
of the accreted/merged gas with any pre-existing gas in the galaxy. 
In either case, we do not know whether the AGN is moving with 
the velocity seen in the ionised gas, or is consistent with the stars. 

We cannot rule out the possibility that 272831 contains two 
SMBHs. In particular, if the observed properties are due to an on¬ 
going merger, if each of the merging galaxies contained an SMBH, 
and if those SMBHs have yet to coalesce, then the two SMBHs will 
currently be present. In this case, the SMBHs may be in a binary 
system or may not yet be under each other’s direct gravitational in¬ 
fluence. However, the assumptions that underpin this scenario can¬ 
not be verified with the available observations. 

A galaxy with similar properties, SDSS J171544.0+600835 
(hereafter J1715), was recently described by Villforth & Hamann 
(2015). In the case of J1715, double-peaked narrow emission lines 
were observed, suggesting a binary SMBH system in which both 
SMBHs are active, and this interpretation was previously confirmed 
by the detection of two distinct X-ray sources (Comerford et al. 
2011). The kinematics of J1715 appear very similar to those of 
272831 (fig. 7 of Villforth & Hamann 2015), in that the stellar and 
ionised gas kinematics are markedly different to each other, with an 
offset of ~50kms -1 in the centre. However, Villforth & Hamann 
(2015) note that the [O III] peaks do not coincide spatially with the 
X-ray sources, indicating that they may not be directly due to the 
two AGN, and suggest a bi-conical outflow as a likely alternative. 
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A counter-rotating gas disc is also suggested, but disfavoured due 
to the ubiquity of outflows in AGN. 

There are some important differences between 272831 and 
J1715: 272831 has much lower [O III] equivalent width and 
[O m]/H/3 ratio, and the rotation curve of 272831 flattens at large 
radii. These differences favour a rotating disc description for the 
ionised gas in 272831, unlike J1715. 

In summary, the IFS observations of 272831 do not provide 
any direct evidence for the existence of a binary SMBH system. 
Under one interpretation of the observations - an ongoing merger 
of two galaxies, each providing one SMBH - we can infer the pres¬ 
ence of a second SMBH in the galaxy, although doing so rests on 
a number of unproven assumptions. In any case, the unusual kine¬ 
matic signature that brought attention to this galaxy is not the result 
of a binary SMBH system, but is related to large-scale motions of 
the gas across the entire galaxy. 

4.2 9016800393 

The bulk of the ionised gas in 9016800393 shows rotational motion 
in the same sense as the stars, with a difference in the maximum ve¬ 
locities due to asymmetric drift. The exception to this ordered mo¬ 
tion is in the centre, where two kinematic components are seen. The 
first component is narrow and has velocity consistent with that of 
the rotating gas and the stars. The emission-line flux ratios indicate 
that this component arises from the NLR of the AGN. However, the 
observed blueshift of the [O III] lines relative to the others indicates 
a more complex kinematic and ionisation structure than expected in 
a classic NLR. 

The second component is broader (304kms _1 ) and 
blueshifted by 230kms _1 along the line of sight. The [Nll]/Ha 
ratio is similar to that of the NLR, but the [Om]/H/3 ratio is 
higher, and the [S ll]/Ha ratio is low. The broadened, blue-shifted 
gas suggests we are observing an outflow directed along or close 
to the line of sight. We note that the relatively strong [O III] 
emission rules out an AGN broad line region (BLR) as the source 
of the broader component, which in any case is significantly 
narrower than the > 1000 km s -1 dispersion typical of BLRs. The 
slight spatial extension of the broader component also favours a 
kpc-scale outflow rather than a BLR or kinematic feature close to 
the accretion disc. 

It is likely that the blueshift of the narrow [O III] lines is a fur¬ 
ther consequence of the complex outflow structure. Comparison of 
the OSSY and two-component LZIFU measurements of the ionised 
gas velocity suggest the OSSY measurement was influenced by the 
strong [O III] lines. Hence, it appears that the identification of an 
offset between the stars and ionised gas was due to the outflow, 
rather than the true velocity of the AGN NLR. 

AGN outflows are known to be ubiquitous, with a range of 
observed properties, and have been shown to produce velocity off¬ 
sets like those seen here (Veilleux et al. 1994, 2005). For example, 
McElroy et al. (2015) found complex emission features indicating 
winds in each of the 17 luminous Type-2 AGN they observed. Of 
their AGN, J103600, J103915, and J152637 are particularly close 
to 9016800393 in terms of their line ratios and velocity dispersions. 

The SFR inferred from stellar populations in 9016800393 
shows a recent burst, with a peak and subsequent decline within the 
last 10 8 yrs. Comparison to the Ha SFR suggests this decline has 
continued: in the outer regions of the galaxy, the current SFR mea¬ 
sured from Ha, 0.060 ± 0.007 M 0 yr _1 , is a factor of ~10 lower 
than the SFR inferred from the youngest stellar populations. We 
speculate that the decline in SFR could be caused by negative feed¬ 


back from the AGN outflow, but testing this scenario requires a 
larger sample of AGN and non-AGN. 

The metallicity of the ionised gas in 9016800393 is typical for 
the galaxy’s stellar mass, so unlike 272831 there is no evidence of 
an external origin for the gas. 

Overall, the observations of 9016800393 strongly disfavour 
a binary SMBH interpretation, and instead point towards an out¬ 
flow from a normal AGN. As noted by CG14 and others, such 
systems are an inevitable contaminant in samples of candidate bi¬ 
nary SMBHs selected by kinematic offsets; spatially resolved spec¬ 
troscopy provides a powerful test for these candidates. 


5 CONCLUSIONS 

We have examined in detail two galaxies that appeared to show 
kinematically offset AGN, using IFS data from the SAMI Galaxy 
Survey. The rich IFS data allow stronger conclusions to be drawn 
regarding the nature of these galaxies. In each case the IFS obser¬ 
vations do not provide any direct evidence for the existence of a 
binary SMBH system. 

In the first galaxy, 272831, there is a strong discrepancy be¬ 
tween the kinematic properties of the ionised gas and those of the 
stars, suggesting a recent merger or gas accretion event. The stel¬ 
lar populations suggest that the star formation associated with the 
event began only 0.1 Gyr ago. The ionisation properties in the cen¬ 
tre of the galaxy are consistent with shock heating, which could 
be due to the AGN or simply the result of accreting gas onto the 
galaxy. 

In the second galaxy, 9016800393, the AGN itself is kinemat¬ 
ically consistent with the host galaxy. However, it is powering an 
outflow that appears blueshifted along the line of sight to the centre 
of the galaxy, producing the observed kinematic offset. 

In both cases, we find no direct evidence for the existence of 
a second SMBH. In the first case we cannot rule out the possibility 
of a second SMBH, but the observed kinematic offset is not due to 
a binary system, and inferring the existence of the second SMBH 
would rely on a number of unproven assumptions. 
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